Crystalline Order on Catenoidal Capillary Bridges 
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We study the defect structure of crystalline particle arrays on negative Gaussian curvature capil- 
lary bridges with vanishing mean curvature (catenoids). The threshold aspect ratio for the appear- 
ance of isolated disclinations is found and the optimal positions for dislocations determined. We 
also discuss the transition from isolated disclinations to scars as particle number and aspect ratio 
are varied. 

PACS numbers: 61.72.Lk Linear defects; dislocations, disclinations 
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Two-dimensional ordered phases of matter on spatially 
curved surfaces have several features not found in the 
corresponding phase for planar or flat space systems [l| . 
For crystalline order on surfaces of spherical topology 
where disclination defects are required by the topology it- 
self, Gaussian curvature can drive the sprouting of discli- 
nation defects from point-like structures to linear grain 
boundary scars which freely terminate in the crystal 

. Even for surfaces such as the torus which admit com- 
pletely defect-free crystalline lattices, the energetics in 
the presence of Gaussian curvature can favor the ap- 
pearance of isolated disclination defects in the ground 
state [1, 0] . For the axisymmetric torus with aspect ra- 
tio between 4 and 10, isolated 5-fold disclinations appear 
near the line of maximal positive Gaussian curvature on 
the outside and isolated 7-fold disclinations appear near 
the line of maximal negative Gaussian curvature on the 
inside [s']. The ground states in these systems are thus 
distinguished by a defect structure that would be energet- 
ically prohibitive in flat space. It is certainly worthwhile 
to explore as many settings as possible in which there are 
qualitative changes in the fundamental structure of the 
ground state, within a given class of order, purely as a 
consequence of spatial curvature. 

The richest confluence of theoretical and experimen- 
tal ideas in the area of curved two-dimensional phases of 
matter has been in colloidal emulsion physics in which 
colloidal particles self-organize at the interface of two 
distinct liquids, either in particle-stabilized (Pickering) 
emulsions Tu] or charge-stablized emulsions [111, [S] ■ 
Two-dimensional (thin-shell) spherical crystals form at 
the surface of droplets held almost perfectly round by 
surface tension. The ordered configurations of particles 
may be imaged with confocal microscopy and the parti- 
cles manipulated with optical tweezers (5|,|3,[l4]. Macro- 
scopic examples of crystalline order on variable positive 
Gaussian curvature surfaces have been constructed by 
forming a soap bubble raft on a spinning liquid [l5| 
and the nature of the order has been analyzed theoreti- 
cally 16]. 

Glassy liquids on negative Gaussian curvature mani- 
folds have also received considerable attention jl7l - l2]| . 
The simplest such manifold conceptually is the constant 
(negative) curvature hyperbolic plane and it even 
appears that the hyperbolic plane can be isometrically 



embedded as a complete subset of Euclidean 3-space, 
although not differentiably [g^. Physical realizations 
of negative Gaussian curvature manifolds in condensed 
matter physics will almost always have variable Gaus- 
sian curvature. The inner wall of the axisymmetric torus 
{S^ X S^) has integrated Gaussian curvature equal to 
—An, balancing an equal and opposite integrated Guas- 
sian curvature on the outer wall. This is responsible for 
the novel ground states noted above. Gaussian bumps 
have regions of both positive and negative Gaussian cur- 
vature and the minimal-type surfaces found in bicontin- 
uous phases of amphiphilic bilayers have spatially ex- 
tended variable Gaussian curvature that is negative on 
average [23j . 

Recently crystalline particle arrays on variable Gaus- 
sian curvature surfaces has been studied by assembling 
particles on capillary bridges formed by glycerol in bulk 
oil spanning two flat parallel plates 24]. Configurations 
may be imaged by confocal microscopy and even ma- 
nipulated with laser tweezers. The interface between the 
inner fluid of the capillary bridge and the outer bulk fluid 
is a surface of revolution with a constant mean curvature 
(CMC) determined by the pressure difference between 
the two fluids [l^, [Igj. Capillary bridges minimize the 
surface area at fixed volume and perimeter and appear in 
the classical work of Delaunay (23, • The value of the 
mean curvature and hence the underlying surface may be 
changed by varying the spacing between the plates. 

Three classes of the Delaunay surfaces have nega- 
tive Gaussian curvature - the nodoids, unduloids and 
catenoids. We will focus on the most analytically 
tractable case of a catenoidal capillary bridge in which 
the mean curvature is everywhere zero. Capillary bridges 
themselves have wide-ranging application. They play 
an essential role in adhesion, antifoaming, the repelling 
coffee-ring effect and in the origin of attractive hydropho- 
bic forces 

The shape of a capillary bridge with mean curvature 
H follows by solving 
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where 2H = ^ + r — r(z) is the representation of a 
surface of revolution with symmetry axis z and Ri and 
R2 are the two principal radii of curvature at any point. 
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The solution corresponding to the special case _ff = (a 
minimal surface) is 

(c cosh(|) cosu \ 
c cosh(^) sin-u , (2) 
V J 

where u £ [0,27r) and v G (— This surface is 
the well-known catenoid parameterized by the radius of 
the waist c located in the z = plane (see FiglT|). From 
the non-zero metric components guu = c^cosh^(^) and 
gvv — cosh^(^), one can obtain the Gaussian curvature 
K = — — ^sech^(-) = — i, where g is the determi- 

Ri R2 ^ c ' 5 ' ^ 

nant of the metric tensor. We see explicitly that the met- 
ric completely determines the Gaussian curvature. The 
solution to Eq.dD) for 7^ is [33] 

/ 7A(6i,t)cosM \ 

x{u,t) = 7A(6l,t)sinu , (3) 

\a + -i{E{e,t) + F{e,t)cose) ) 

where t € [^0,^1], u € [0,27r), E{e,t) = 
Jl;Ai9,i)di, F{e,t) = 1/A{0, i)di and A{e,t) = 

— sin^ 6 sin^ t. 7 plays the role of a scale factor. The 
curves generated are periodic in t with period tt and have 
maxima at t = kir and minima at t = {k + 1/2)'k for in- 
teger k. The value of 6 controls the shape of the profile: 
it generates an unduloid for 9 E [0, 7r/2), a nodoid for 
9 € (7r/2,7r); and a semi-circle for 9 = tt. The shapes of 
the capillary bridges in the work of ^3] can be fit by 
Eq.Q, the general CMC expression with H 0. 

Here we study crystalline order on the simplest case of 
a catenoidal capillary bridge (H = 0) in the framework of 
continuum elasticity theory [l], [3, Hj] • For simplicity, we 
measure all lengths in units of the radius of the contact 
disk. 

The topology of the capillary surfaces we study is that 
of the annulus, with Euler characteristic zero, since the 
liquid bridge makes contact with the plates at the top 
and bottom. Such a surface admits regular triangula- 
tions with all particles having coordination number 6. 
Although defects (non 6-fold coordinated particles) are 
not topologically required they may be preferred in the 
crystalline ground state for purely energetic reasons since 
negative Gaussian curvature will favor the appearance 
of 7-coordinated particles (-1 disclinations) . To deter- 
mine the preferred defect configuration we map the mi- 
croscopic interacting particle problem to the problem of 
discrete interacting defects in a continuum elastic back- 
ground. The defect free energy Fei, in the limit of van- 
ishing core energies, may be expressed in the form p].[l6j 

Fei^l-Y f G2L{x,y)p{x)p{y) d'xd^y. (4) 

Here G2L {x, y) is the Green's function for the covariant 
biharmonic operator on the surface M, Y is the Young's 
modulus for the crystalline packing, and x and y are po- 
sition vectors on the surface. The effective topological 
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FIG. 1: (a) The three dimensional shape of a catenoid with 
aspect ratio c = 0.85. (b) A catenoid with aspect ratio c — 0.9 
(green) deforms to c = 0.7 (blue). 

charge density is p{x) = ^q{x) — j^^Wjbi(x) — K{x), in 
which the first term is the disclination charge density 
lix) = 9c(^(^ ~ Xa), the second term is the dislo- 
cation density bi{x) = X]^^f'^('^~ ^p) ^^'^ K{x) is the 
Gaussian curvature. In this expression 7'-' — e^^ 
By introducing x ^^^d F such that A^x(a;) = Y p{x) and 
r(a;) = Ax(a;), Eq.Q can be written in a more compact 
form 

Fei = ^ j^T\x)dx, (5) 

in which V(x)/Y — j GL(x,y)p{y) dy + U(x) with 
l\U{x) = and GL{x,y) is the Green's function for the 
covariant Laplacian on M which satisfies l\GL{x,y) — 
S{x, y)]x,y £ M with the boundary condition Gl{x, y) — 
0',x,y € Dm- By conformally mapping the surface param- 
eterized by {u, v} onto an annulus in the complex plane 
via z = p{u, u)e™, the Green's function Gl is found to be 
M Gi(x,y) = 2Lin|(p-iz(f)_z(y))/(l-p-iz(f)z(y))|, 
in which po is the radius of the outer circle of the annu- 
lus in the complex plane. For a catenoid, the confor- 
mal mapping is given by p{u,v) — c e'"'/^ and pa = 
c exp(arcsech(c)). 

Disclinations are expected to appear in the crystalline 
ground state when the Gaussian curvature is sufiicient 
to support them. Consider therefore a putative isolated 
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FIG. 2: (a) A family of geodesies in {u, v} coordinates cen- 
tered at a point on the waist of a catenoid with c = 0.85. 
(b) The \K{0)\-!vr^ (dashed curve) and numerical result (solid 
curve) of the integrated Gaussian curvature over a geodesic 
disk of radius r versus r. c — 1/2. 



disclination of strength q = —1 (coordination number 
7) at the waist of a catenoid. The curvature condition 
above requires that there exist a disk of geodesic radius 
Tc, centered on the 7-disclinaton, for which [H, 113 



KdA = . 

disk 3 



(6) 



Clearly Vc must be less than the geodesic distance I from 
the waist to the boundary . For a given size catenoid 
c, we calculate I and the integral of the Gaussian curva- 
ture over the geodesic disk of radius I. The value of c for 
which the integrated curvature equals —n/3 is the criti- 
cal value of c for the appearance of 7-disclinations. We 
compute the integral of the Gaussian curvature numeri- 
cally. We first construct a family of geodesies radiating 
from the core 7-disclination {a.t u — 0,v — 0) by solving 
the geodesic equation: 
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dxP dx' _ 
d\ dX 



and 
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(7) 



is the Christoffel 



in which x = u,x = v 

symbol of the second kind. This second order differen- 
tial equation has a unique solution given an initial po- 
sition and an initial velocity. The initial conditions are 
xi(0) = a;2(0) = 0, i|i|o = IcosO, and ^|o = sm0, 
where is the angle of the initial velocity with respect 
to e". Given a geodesic radius r, the coordinates of the 
end point of the geodesic curve can be found. These end 
points form the boundary of a disk in {u, v} coordinates 
(see FigUJa)). We then integrate the Gaussian curva- 
ture over the prescribed disk numerically. The critical 



value of c is found to be c* = 0.85 and the corresponding 
critical radius is Tc = 0.53. Note that integrated Gaus- 
sian curvature for this critical catenoid is quite large [l^l : 
J KdA ~ —6.6. The critical value c* can also be esti- 
mated as follows. By introducing Gaussian normal co- 
ordinates (r, 9) centered on a 7-disclination at height zq 
above or below the waist of the catenoid, the effective 
(screened) disclination charge at distance r is given by 

mil 
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7r_sech^( — ) + 0{r'). 
d- c 

The critical radius is reached when the effective disclina- 
tion density vanishes: Pefjij'c) — 0. For a 7-disclnation 
on the waist (zq = 0) this gives r^jc = dc~ ■\/l/3 w 33°. 
Now on the catenoid the geodesic length from the waist to 
the boundary is " cosh(v/c)dw — Vl — c^- The critical 

catenoid size c* is then given by Tc* Vl — c*^. This 
yields c* = ■\/3/2 w 0.87. This estimate for c* is very 
close to the numerical value 0.85. Why are these two 
values so close? In calculating the effective disclination 
charge, we use Kifi)'r:r'^ to approximate the integral of 
the Gaussian curvature over a geodesic disk of radius r. 
The Gaussian curvature is overestimated as its magni- 
tude is maximum at r = (on the waist). On the other 
hand, since K{Q) = lirrir^o^^i'^T'^ — A{r)) < 0, the real 
area A{r) of the disk with geodesic radius r is bigger 
than 7rr^, i.e., the disk area is underestimated in our ap- 
proximation. These two approximations tend to cancel 
each other out. For a typical value of c = 1/2, |i<r(0)|7rr^ 
and the numerical result of the integral of the Gaussian 
curvature versus r is plotted in Figj^l^b). As expected 
the flat space approximation K{0)Trr'^ is good for small 
r(r < 0.2). 

The critical waist size c* can also be estimated from 
energetic arguments. From the free energy of Eq.([5|) one 
can analyze the geometric potential describing the inter- 
action between disclinations and the intrinsic Gaussian 
curvature of the surface. The result is shown in Figl3|). 
We see that the optimal position of a disclination shifts 
from the boundary to the waist as c decreases. The tran- 
sition point for the emergence of a disclination in the 
interior of a catenoidal capillary bridge is c* « 0.8, again 
consistent with the value obtained above based on geo- 
metrical arguments. 

Net disclination charges may appear either in the 
form of point-like isolated disclinations or extended lin- 
ear grain boundary scars. Scars result from the screening 
of an isolated disclination by chains of dislocations and 
typically arise when the number of particles exceeds a 
threshold value beyond which the energy gained exceeds 
the cost of creating excess defects Here we semi- 
quantitatively construct the phase diagram for isolated 
disclinations versus scars on a catenoidal capillary bridge 
characterized by the number of particles and the aspect 
ratio of the catenoid c. 
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FIG. 3: (a) The geometric potential of an isolated disclination 
for three different values of c: c = 0.8 (red), c = 0.75 (green) 
and c = 0.7 (blue). The optimal position of an isolated discli- 
nation moves from the boundary to the waist of the catenoid 
in the rather narrow window c between 0.8 and 0.75. (b) The 
geometric potential of an isolated disclination for catenoids 
with c = 0.5 (red) and c = 0.6 (green). 
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FIG. 5: (a) The geometric potential of isolated dislocations 
as a function of height for four different values of c: c = 0.7 
(black), c — 0.68 (red), c = 0.66 (blue) and c = 0.65 (green). 
The optimal position moves from near the boundary towards 
the waist as c decreases, (b) is an enlarged view of the blue 
curve to show the transition from the red to the green curve. 
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FIG. 4: The phase diagram in the particle number-aspect 
ratio plane for isolated disclinations versus scars for c < c*. 



Consider a disclination on a capillary bridge (for c < 
0.85) radiating m grain boundaries (scars). The spacing 
of neighboring dislocations is Z = a m/seff [2j, where a 
is the lattice spacing. As s^fj — > 0, the dislocation spac- 
ing within a scar diverges and the grain boundary termi- 
nates. If the disclination can be completely screened by 
Gaussian curvature within a circle of radius r « 3a, then 
grain boundaries will not form around the core disclina- 
tion. The condition for isolated disclinations is therefore 
\KmaxTTiSa)'^\ 7r/3, whcrc \Kjnax\ = is the Gaus- 
sian curvature at the waist of the bridge. On the other 
hand, the number of particles N is related to the sur- 
face area A between z S [—Zm,Zm] via A{c) = ^a^N. 



The curve separating isolated disclinations from scars is 
thus given by = i8v^^A(c) ^ plotted in Figgja). The 
phase boundary reveals two basic types of transition in 
the topological structure of the ground state as the par- 
ticle number and the geometry (aspect ratio) of the cap- 
illary bridge are varies. For a fixed catenoid aspect ratio 
below the critical value for the appearance of excess 7s in 
the interior there is a transition from isolated 7s to linear 
grain boundary scars with one excess 7 as the number of 
particles increases. For a fixed number of particles above 
a threshold value {Nc ~ 300) there is a transition from 
isolated disclinations to scars as the capillary bridge gets 
fatter and the decreasing Gaussian curvature is insuffi- 
cient to support isolated 7-disclinations. 

Disclinations and anti-disclinations attract and may 
form dipole bound states (7-5 pairs). Such dipole config- 
urations are themselves another type of point-like topo- 
logical defect in two-dimensional crystals - dislocations. 
Dislocations on a triangular lattice correspond to two 
semi-infinite Bragg rows 60° apart both terminating at 
a common point - the location of the dislocation. Since 
they are tightly bound states of disclinations the ener- 
getics of dislocations may be derived from the govern- 
ing energetics of disclinations on a curved geometry - 
Eq.([5]). Dislocations, unlike disclinations, are oriented. 
The Burgers vector characterizing a dislocation at po- 
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FIG. 6: The interaction of defects, (a) The dislocation- 
dislocation interaction Vdd of two dislocations along e" at 
the same height, z = 0.5 (black), z — 0.3 Zm (red) and 
z = 0.1 Zm (green), (b) The disclination-dislocation interac- 
tion Vod as a function of their longitudinal separation. The 
disclination is fixed on the waist (c = 1/2). 

sition x" is perpendicular to the 5-7 bond. An analysis of 
Eq. ([5]) shows that the preferred orientation of the Burg- 
ers vector is along e". This is clear from the fact that 
the 7-disclination has minimum energy when located at 
the waist with the accompanying 5-disclination in the 
direction of the boundary where the negative Gaussian 
curvature drops most rapidly. Thus the 7-5 bond should 
along a meridian and the Burgers vector along a line of 
latitude. From here on we restrict ourselves to this case. 

The variable Gaussian curvature of a catenoidal capil- 
lary bridge also leads to optimal positions for isolated 
dislocations. FiglS] shows the geometric potential for 
isolated dislocations as a function of height above the 
waist. As the waist radius c decreases the optimal po- 
sition of an isolated dislocation moves from the bound- 
ary to the interior of the capillary bridge since the in- 
creasing maximal negative Gaussian curvature increas- 



ingly attracts 7-disclinations with their tightly bound 
5-disclinations. The boundary-to-interior transition oc- 
curs for c** w 0.68. The corresponding integrated Gaus- 
sian curvature J KdA « —9.0 fl^l- This detachment 
transition is also observed in experiments with capil- 
lary bridges - in the experimental case the capillary 
bridges are generally nodoids with non-vanishing mean 
curvature and the analysis is corresponding more elab- 
orate. The optimal position of a single dislocation for 
small c, say c = 0.3, is z{c — 0.3) — 0.25zm- Thus 
the optimal position of a single dislocation is far from 
the waist, even for a strongly curved catenoidal capillary 
bridge, in contrast to the case of disclinations. This result 
can be understood in terms of the Peach-Koehler forces 
acting on the individual positive and negative disclina- 
tions that make up a dislocation [23, • While the 7- 
disclination prefers to be at the waist the 5 prefers to be 
at the boundary - the competition results in an optimal 
dislocation position somewhere in between the two ex- 
tremes. Here we treat only single dislocations but it is 
possible for chains of dislocations to appear in the form 
of pleats, as elegantly discussed in Ref . j24| . 

Finally we turn to the interaction between defects 
themselves. FiglH] (a) shows the dislocation-dislocation 
interaction along e". Two dislocations at the same height 
feel a short-range repulsion. Note that near the waist 
shallow local minima appear. This differs from the in- 
teraction in fiat space where parallel dislocations always 
repel to each other with a logarithmic potential [37|. 
The attractive interaction between a disclination and a 
nearby dislocation is shown in FiglS^b) with the disclina- 
tion fixed on the waist. 

The influence of spatial curvature and topology on 
two-dimensional phases of matter continues to yield sur- 
prises. The presence of 7-disclinations in negative curva- 
ture crystals may offer unique opportunities for function- 
alization of micron-scale crystallized "superatoms" via 
chemistry that recognizes the unique crowded environ- 
ment of a 7-disclination [l], [ll] . 
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